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The near-forward cross section for the pp — ► pnir + reaction has been measured 
at 492 MeV by using the large acceptance ANKE magnetic spectrometer placed 
at an internal target position of the storage ring COSY-Julich. Protons and pions 
emitted near zero degrees were detected in coincidence, and those with 9^ < 2° and 
9 P < 2.5° were subjected to detailed analysis. Under these conditions, the excitation 
energy in the np system was below 3 MeV over the measured momentum range. 
This is the region of the np final-state-interaction peak, which was measured with 
a resolution of a fraction of an MeV. The shape of the peak allows one to conclude 
that the fraction of final spin-singlet np pairs is below about 10%. By using the 
results of scattering theory, this limit is confirmed through a comparison with the 
cross section for pp — > dir + . The smallness of the singlet contribution is consistent 
with trends seen in lower energy data. 
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The two-body pp — > dir + reaction has long been used as a test bed for both ex- 
perimental measurements and theoretical models of intermediate energy pion 
production [1]. Unfortunately, there is no 1 Sq bound state in the S-wave spin-0 
np system, analogous to the 3 Si deuteron bound state in the spin-1 system, 
that would enable a similar study of the other isospin channel. On the other 
hand, there is a virtual (antibound) 1 5'o state very close to threshold at an 
np excitation energy E np xs —0.07 MeV. Both this and the deuteron pole 
at —2.22 MeV induce strong S'-wave final-state-interaction (fsi) peaks in the 
pp — > pnn + reaction leading to spin-singlet and triplet np channels. To sim- 
ulate a two-body reaction, the excitation energy in the np system should be 
confined to within a few MeV. Furthermore, the singlet and triplet contribu- 
tions should be cleanly separated. This should in fact be possible with good 
energy resolution because the fsi peak widths are proportional to the (real or 
virtual) binding energies. The width of the singlet peak would therefore be a 
fraction of an MeV as compared to the 2-3 MeV of the triplet. 

Inclusive pp — > n + X reactions at intermediate energies show an np enhance- 
ment at the edge of phase space, but the energy resolution achieved and the 
contamination from the much stronger pp — > dii + channel makes it very diffi- 
cult to study the fsi region at low E np [2]. The data nevertheless suggest that 
the triplet fsi provides the major signal [3]. To avoid feed-through from the 
dn + final state, it is useful to measure the proton and pion of the pnn + channel 
in coincidence and this can also lead to better determination of E np . Such an 
experiment was carried out at LAMPF at 800 MeV [4]. Although fsi peaks 
were observed, the number of E np points in the very low energy region was 
limited, as was the resolution, and so it was not possible to isolate a narrow 
singlet fsi peak from a broader triplet one. Nevertheless their model-dependent 
analysis suggested that the singlet fsi represented a very large fraction of the 
total strength at all angles measured. We report here on the results of an 
exclusive pp — > pnir + measurement at T p = 492 MeV, where both the proton 
and pion were detected close to the forward direction. The resolution of a 
fraction of an MeV in E np allows us to separate contributions from the singlet 
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and triplet fsi peaks and our results suggest that the singlet fraction is below 
10%. This upper bound is confirmed by comparing the magnitude of our cross 
section with that for pp — > dn + using general scattering theory results [5]. 

The experiment was carried out at the ANKE spectrometer placed inside 
the storage ring of the proton synchrotron COSY-Jiilich. The main ANKE 
components are shown in Fig. 1, but see Ref. [6] for a detailed description. 
Foil targets made of CH2 (polyethylene) were used in the experiment and 
the background from carbon was measured using a polycrystalline diamond 
target. The spectrometer dipole D2 separates positively charged ejectiles from 
the circulating COSY beam and allows one to determine their momenta at 
the target over a wide momentum range. The horizontal angular acceptance 
of |i?h| < 12° is mainly fixed by the allowed start-stop combinations of the 
on-line TOF-trigger system. The vertical angular acceptance is in the range 
|$v| < 7° (3.5°) for the lowest (highest) detectable ejectile momenta. 
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Fig. 1. Top view of the ANKE spectrometer and detectors used for the FSI studies. 

The detection system [7] comprises two multi-wire proportional chambers 
(MWPC 1, 2 in Fig. 1) for track and momentum reconstruction, 23 scin- 
tillation counters close to the vacuum window of D2 for time-of-flight start 
measurements and 15 range telescopes located along the focal surface of D2. 
The side wall counters allow larger ejectile momenta to be measured. The 
identification of pions and protons is achieved via TOF and energy loss mea- 



3 



surements. Background originating from the pole shoes of D2 is suppressed 
with the help of the track information from the MWPC's. The field strength 
of B = 0.861 T in D2 was chosen so that the pions used in the analysis 
were detected and identified in telescope #14 and the protons in the side- wall 
counters. 

The spectrum of missing-mass-squared for the p(p,pir + )X reaction on a CH 2 
target, shown in Fig. 2 for the full ANKE angular acceptance, displays a 
prominent neutron peak, with a standard deviation of cr{m^) = 0.0034 m^, 
which is not present for the C target. The integrated luminosities for the two 
targets are different but, when the relative CH 2 /C normalisation is evaluated 
from the spectra at high missing masses, it is found that the neutron peak 
sits on a carbon background of a few percent, which can easily be subtracted. 
Fits to the position of the deuteron and neutron peaks determined the proton 
beam energy to be T p = (492 ±1.7) MeV. 
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Fig. 2. Spectrum of the square of the missing mass for the p(p,p7r + )X reaction, in 
units of the neutron mass, for: a) CH2; b) C; c) CH2— C targets, normalised to the 
same luminosity. By comparing the spectra at high masses, the relative luminosity 
L(C)/L(CH 2 ) is deduced to be 0.264 ± 0.020. 

In order to obtain clean samples of events with low np excitation energies, 
software cuts were imposed so as to extract two groups of events. Both of 
these have pion angles 6^ < 2.0°, but the first group includes protons with 
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angles 6 P < 2.0°, whereas the proton cut was extended to 2.5° for the second 
set. After subtracting the C background, the number of events remaining in 
the two cases was 1312 and 2008 respectively. The ratio follows that of the 
proton solid angle, the cut in the pion angle being of no importance here. 

In the raw proton momentum spectra shown in Fig. 3, the carbon background 
with our cuts is below 3%. The central np excitation energies are indicated 
and, from this scale, it is seen that our events generally correspond to E np < 
3 MeV, which means that the fsi region is well covered. Because of the limited 
statistics, we have summed all events over the whole angular ranges and this 
introduces an uncertainty of AE np 0.55^ E np MeV. This is far greater than 
the intrinsic resolution of the system which, through kinematic fits to the 
pp — > pnn + reaction, is about a = 160 keV. Thus the width of any spin- 
singlet contribution in Fig. 3 is determined by the angular integration rather 
than the natural width or the intrinsic energy resolution of the apparatus. 

With the angular cuts we have imposed, the proton measurement efficiency is 
about 60% at small proton momenta (< 350 MeV/c), constant at about 99% 
for momenta between 350 MeV/c and 410 MeV/c, but falls rapidly thereafter. 
Taking this and other effects into account in a detailed Monte Carlo study, 
we show in Fig. 3 a simulation where we have assumed that the cross section 
is proportional to three-body phase space times essentially a pure spin-triplet 
fsi factor. The significance of the excellent agreement with the shape of the 
spectrum will be crucial for the later discussion. One of the advantages of our 
experimental set-up is that we also have a measurement of the correspond- 
ing pion momentum spectrum. However, within the current precisions, the 
singlet /triplet ratio is better determined from the proton spectrum [8]. 

The absolute normalisation of the pp — * pnii + cross section is achieved by 
comparing the pp — > dn + events, measured in parallel, to standard cross sec- 
tion compilations [1]. Because of the dominance of the dir + final state, the 
reaction could be identified with sufficient precision without having to detect 
the deuteron, thus avoiding uncertainties arising from deuteron break-up in 
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Fig. 3. Raw proton-momentum spectra of pp — > pn7r + events (points with statistical 
errors), for events with 0^ < 2°. All the histograms are normalised to the total 
number of events. The relative np energy is indicated for protons and pions emitted 
at 0°. Figure a) is for 9 P < 2° and b) for 6 P < 2.5°. The dashed histogram is a 
Monte Carlo simulation assuming pure spin-triplet final state (£ = 0). Fits to the 
fsi formula of eq. (2) are consistent with this, giving £ = —.008 ± .038 for case a) 
and £ = —.008 ± .036 for b. A statistical mixture of spin states, i.e. £ = 0.25, leads 
to the dotted histograms with x 2 /ndf =151.2/17 for a) and x 2 / nd f =178.2/17 for 
b). 



the counters etc. The contamination from the carbon background is here larger 
(~ 23%) than for three-body events, but it is easy to correct for this using 
the normalised results from the carbon target discussed earlier. Our results 
are presented in Fig. 4, where only statistical errors are shown. The principal 
systematic effects are: 

• Uncertainty in the number of pp — > dn + events, due to the tail of pions from 
the pp — > pnir + reaction where the proton escapes detection: pa 7%. 

• Uncertainty in the ratio of C/CH 2 integrated luminosities: w 2.5%. 

• Uncertainty in the effective proton solid angle: ps 1.5%. 
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• Uncertainty in the SAID solution for pp — > dn + , used for normalisation: 
« 2.0%. 

Summing these uncertainties quadratically, we arrive at an estimate of the 
total systematic error of 8%. 




Fig. 4. Five-fold differential cross section for the pp — > pn7r + reaction as a function of 
the measured proton momentum for events with 6 n < 2°. The np excitation energy 
is indicated for protons and pions emitted at 0°. The upper figure is for 6 p < 2° and 
the lower for B v < 2.5°. The solid curves correspond to predictions for the spin-triplet 
contribution to the fsi peak, obtained from eq. (3) using pp — > dir + input data from 
the SAID SP96 solution [1], after averaging over the angular acceptance. The dashed 
curves are raw predictions without averaging. The estimates are expected to be only 
weakly model-dependent [5]. 



Since the pp — > pmr + events of Fig. 3 and Fig. 4 are all concentrated in the 
region E np < 3 MeV, we can make an estimate of the relative amounts of 
spin-singlet and -triplet final states using the fact that any singlet fsi peak 
should be so much narrower than the triplet. We assume a variant of the 
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Goldberger- Watson fsi formula [9] 

FSM*) = Aig±g, (1) 

where the ctj and fa are determined from low energy np scattering data [13]. 
The triplet and singlet parameters are a t = 0.232 fm -1 , a s = —0.040 fm -1 , 
= 0.91 fm -1 , and (3 S = 0.79 fm -1 . It should be noted that the numerators in 
eq. (1) vary very little over our restricted E np range, so that the fits are weakly 
dependent upon the values of the The overall constant N { is determined 
by the condition that the integral of FSIj(A;) over k for E np < 3 MeV is unity, 
which requires that N t = 0.104 fm and N s = 0.0297 fm. We then fit our data 
with the form 

d5 ^ (P ^^ 7r+) oc (Kinematic factor) x [£FSI s (Jfe) + (1 - FSI t (Jfe)] , 

KXjJ'P C_li Lip Lli Lq^ 

(2) 

so that £ has the physical significance of being the fraction of the cross section 
leading to singlet np states for E np < 3 MeV. 

When this form is passed through the Monte Carlo simulation and normalised 
to the data in Fig. 3, an excellent fit is achieved with £ = —.008 ± .038 and 
X 2 /ndf= 18.2/17 for group a) and £ = -.008 ±.036 and x 2 /ndf= 15.9/17 for 
group b). Effects of smearing the predictions over the spread in E np induced 
by the angular acceptance and momentum bite etc. have been taken into 
account. A value of £ = 0.25, which would correspond to a purely statistical 
factor, gives excess of events in the region around zero excitation energy and 
shortage at the right tail of the spectrum. This increases the x 2 / n df by about 
10 units. Thus we can assert purely from the shape of the fsi peak that the 
singlet contribution to the cross section is likely to be below 10%. 

A second proof that the spin-singlet fraction is small can be found from the 
normalisation rather than the shape of the measured differential cross sec- 
tion. The Faldt- Wilkin extrapolation theorem links the np scattering wave 
function to the deuteron bound-state function, independently of the form of 
the np potential [10,3], and this allows one to predict the S'-wave spin-triplet 
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contribution to the pp — > pmr + cross section at low E np in terms of that for 
pp — > dir + [5]. To the extent that deuteron D-state effects can be neglected, 

" t7^Z + ' = (Kinematic factOT) x 0?T^|) x 5^ "* * +) ■ 

(3) 

This relation does not require explicit knowledge of the pion production op- 
erator, merely that it be of short range. It describes well [5] the fsi region of 
the LAMPF pp — > pn-n + data taken at a variety of angular configurations [4], 
though it should be stressed that the number of experimental points in the fsi 
peaks was small. 

Because eq. (3) has the same dominant fsi factor as that of the spin-triplet 
Goldberger- Watson factor of eq. (2), it is no surprise that the shape of the 
differential cross section as a function of the final proton momentum shown in 
Fig. 4 is largely reproduced. To within possible corrections of the order of the 
deuteron D-state probability (~ 6%), this approach provides a robust lower 
bound on the three-body cross section. Although the curves pass through the 
points to the left of the peak, they are slightly low compared to the data on 
the right. This is compatible with theoretical and experimental uncertainties, 
and so the normalisation of the data also shows that any singlet contribution 
must be below about 10%. 

The inclusive TRIUMF p(p, ir + )X data [2] indicate significant contributions 
from np final states other than triplet S-wave, but their excitation energies are 
such that P-waves cannot be neglected. The exclusive p(p,p7r + )n measurement 
carried out at CELSIUS at 400 MeV [11] was analysed as for a single-arm 
experiment by integrating over the final proton momentum as well as over 
the pion angle. When this is done, there is an indication of contributions 
other than spin-triplet S'-waves, though there is only one high point in the 
E np < 3 MeV range. 

Given that the singlet fraction is so small at intermediate energies, this sug- 
gests that one should measure its production directly by studying the pp — > 
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pp7T° reaction, which could be done in parallel with pmr + detection at ANKE. 
The reaction has already been investigated at CELSIUS up to 425 MeV and 
values of the differential cross section quoted separately for E pp < 2.6 MeV [12]. 
Though there is no similar detailed measurement of the pp — > pnn + cross sec- 
tion with which to compare the CELSIUS tt° data, it is possible to estimate 
the spin-triplet contribution quite reliably from that of pp — > d-n + by employ- 
ing the Faldt- Wilkin theorem [3]. If then the cross section ratio is extrapolated 
to 492 MeV, it indicates that the singlet fraction should be below about 1%. 
This very low value is due, in part, to n + production being maximal in the 
forward direction [1], whereas the 7r° are produced preferentially at 90° [12]. 
Thus the singlet fraction at 492 MeV would be expected to be much larger for 
an experiment carried out at 90°. 

In summary, we have measured the exclusive pp — > pnn + cross section reaction 
near the forward direction with a resolution in E np superior to that of previous 
works. This allows us to put an upper bound on the proportion of singlet np 
final states, which is confirmed independently by comparing with pp — > dir + 
results. It is hoped to continue the investigation at other energies, studying 
simultaneously the pnn + and ppn° channel. 

The results presented here were obtained during the ANKE commissioning 
experiments. We would like to thank the colleagues who helped to realise the 
ANKE spectrometer, in particular O.W.B. Schult, K. Sistemich and R. Maier, 
and the many former Diploma and Ph.D. students. We express our thanks 
to the COSY team for its dedication and very smooth implementation of 
ANKE into the accelerator ring. During the measurements, support from the 
"Zentrallabor fur Elektronik" (ZEL) of the FZJ, in particular W. Erven, was 
indispensable. 
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